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ABSTRACT 

We present a new method for recovering the kinetic temperature of the intervening 
diffuse gas to an accuracy of 10%. The method is based on the comparison of 
unsaturated absorption-line profiles of two species with different atomic weights. 
The species are assumed to have the same temperature and bulk motion within the 
absorbing region. The computational technique involves the Fourier transform of 
the absorption profiles and the consequent Entropy-Regularized ^-Minimization 
[ERM] to estimate the model parameters. The procedure is tested using synthetic 
spectra of C + , Si + and Fe + ions. The comparison with the standard Voigt fitting 
analysis is performed and it is shown that the Voigt deconvolution of the complex 
absorption-line profiles may result in estimated temperatures which are not physical. 
We also successfully analyze Keck telescope spectra of C IIA1334 and Si IIA1260 lines 
observed at the redshift z = 3.572 toward the quasar Q1937-1009 by Tytler et al.. 



Subject headings: cosmology: observations — methods: data analysis — quasars: 
absorption lines — quasars: individual (1937-1009) 



1 Based in part on data obtained at the W. M. Keck Observatory, which is jointly operated by the University of 
California and the California Institute of Technology. 
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1. Introduction 

New generation of large telescopes enables to study the intergalactic and interstellar absorbing 
gas with high instrumental resolutions sufficient to resolve absorption features into the individual 
subcomponents. From absorption-line spectra, one usually estimates the column densities of 
absorbing species, N, the Doppler 6-parameter of the diffuse gas, and the positions (i.e. radial 
velocities) of the subcomponents which are traditionally considered as being caused by the 
intervening 'clouds'. 

The observations show the ubiquitous asymmetry of the line profiles (Spitzer & Fitzpatrick 
1993; Huang et al. 1995; Lu et al. 1996; Prochaska & Wolfe 1997). Their widths are believed to 
reflect the sum of thermal and non thermal ('bulk' or 'turbulent') motion. It should be noted 
that the term 'turbulent' used in spectroscopy is less distinct as compared with hydrodynamic 
turbulence and labels the lack of information about the nature of the broadening mechanism 
(Traving 1975). Therefore to treat observations a model of the line broadening is to be specified 
first of all. 

It is a commonly used assumption that both thermal and turbulent velocities along a given 
line of sight can be represented by Gaussian distributions, i.e. thermal and non thermal broadening 
can be added quadratically 

b 2 = v 2 h + 2a 2 t . (1) 
Here at is the rms turbulent velocity and Vth is the thermal velocity 

v th = ^2fc B T/m a , (2) 
where k-g is Boltzman's constant, m a is the mass of the element, and T the kinetic temperature. 

To separate the contribution of thermal and turbulent motion one usually utilizes the line 
profiles of different elements assuming that they have the same temperature and the velocity 
distribution within the individual components. If this is the case, then higher mass elements 
should have lower apparent 6-values provided the measurements have a sufficiently high spectral 
resolution. With such reservations one can solve equation (|l|) for T and at for different species. 

Applied to the real absorption-line spectra, this technique often leads, however, to the 
implausible wide range of T among individual 'cloudlets' yielding for some of them even unphysical 
negative kinetic temperatures - clearly an absurdity (Spitzer & Fitzpatrick 1993, 1995; Fitzpatrick 
&; Spitzer 1997; Tripp, Lu & Savage 1997). This suggests that in general the problem of the 
6-parameter measurement from the interstellar and/or quasar metal absorption-line spectra is far 
from being solved. 

To a certain extent this problem is connected with the interpretation of complex absorption 
features as being caused mainly by the density fluctuations. However, as shown by Levshakov & 
Kegel (1997) the intensity fluctuations within the line profile can occur due to correlation effects 
in the large scale velocity field even for the homogeneous gas density. 
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The present paper is primarily aimed at the problem of the kinetic temperature measurement 
accounting for these correlation effects. The analysis employs absorption profiles of two species 
with different atomic weights. These profiles are considered to arise in a transparent absorbing 
region of thickness L with a correlated velocity field and a homogeneous temperature T. Although 
in principle density fluctuations can be incorporated, here we confine ourselves to the case of the 
homogeneous gas density and, hence, the density ratio of the chosen elements is constant along 
the line of sight. These conditions yield the similarity in the observed profiles of different elements 
(examples can be found in Spitzer & Fitzpatrick 1993, or in Songaila 1998). The assumption of 
constant T is supported by the fact that the states of thermally stable, optically thin interstellar 
gas take only limited temperature ranges over a relatively wide range of density variations (see, 
e.g., Figure 8.9 in Shu 1992). 

The Entropy-Regularized ^-Minimization (the ERM method) we propose allows a robust 
measurement of a single value of T from the multicomponent (complex) absorption spectra in 
cases where the standard Voigt-profile fitting procedure yields different values of T for the different 
subcomponents. The description of our model and its basic assumptions is given in § 2. A study 
of the accuracy of the method is provided through a series of numerical simulations of complex 
spectra produced in the correlated velocity fields of different structure which is presented in § 3. 
We proceed further in this approach with high resolution C IIA1334 and Fe IIA2600 synthetic 
spectra in § 4, in order to derive the gas temperature and to compare its value with the standard 
Voigt fitting procedure yield. An example based on real observations of high-redshifted metal 
absorption lines is presented in this section as well. We summarize our conclusions in § 5. 



2. Data analysis technique 

In this section we present a compendium of the standard radiative transfer equations 
describing pure absorption in spectral lines, and the consecutive analytical technique to derive the 
T-estimation and the radial- velocity distribution function </> t from the line profiles of two elements 
with different atomic weights. Here, we are dealing with noiseless data, assuming that the given 
spectra are exact. 



2.1. Basic equations 

If the interstellar (intergalactic) absorption line is optically thin, one observes directly the 
apparent optical depth r* (A) as function of wavelength A 

T*(A) = ln[/ (A)// obs (A)] , (3) 
where I bsW and -^o(^) are the intensities with and without the absorption, respectively. 
The recorded spectrum is a convolution of the true spectrum and the spectrograph 
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point-spread function sp 

/+oo 
Jo(A / )e-^ A ')0 sp (A-A')dA / , (4) 
-oo 

where r(A) is the true (intrinsic) optical depth. 

The intensity of the background continuum source Io(\) changes usually very slowly over the 
width of the spectrograph function, and hence (||) becomes 

/+oo 
e- T ( A ')0 sp (A-A / )dA'. (5) 
-oo 

For the following it is convenient to change from A-space to the radial-velocity v-space. 
Further, let the line-absorption coefficient reckoned per atom (the absorption cross-section) at 
radial velocity v within the line profile be denoted by k(v). Then the optical depth at velocity v is 

T (v) = Nk(v). (6) 



For the case of unsaturated absorption lines when the effect of radiation damping is negligible 
k(v) may be described by the Doppler kernel of the Voigt function (Spitzer 1978). Then the 
profile function is defined locally by the thermal broadening, i.e. the linewidth is determined by 
the kinetic temperature only (b = t>th)- Now k(v) can be expressed as the product of a constant 
ko for a particular line and the Gaussian line profile function 4>th( v ) 

k(v) = k ^ th (v) = — / abs A -j^e^-^K , (7) 
m e c Vvru th 

where m e and e are the mass and charge of the electron, / a b s and Ao are the oscillator strength 
and the wavelength of the line center at which the corresponding radial velocity is equal to wo- 
lf atoms are involved in both thermal and hydrodynamic motions and the temperature is 
constant along the line of sight, then the total line profile function is a convolution of (fithiv) an d 
the radial- velocity distribution function (pt(v) caused by bulk motion 

/+oo 
M^th(v-v')dv' , (8) 
-oo 

and t(v) = Nk ^ to t(v). 

Both functions 4>t and 0th are normalized so that J(-)dv = 1. In principle, 0t reflects the 
density fluctuations as well, and in reality is the density weighted radial- velocity distribution. But 
while we are dealing with the similar profiles the explicit treatment of the density fluctuations can 
be avoided. 

The instrumental point-spread function sp can be determined experimentally, whereas <pt 
and ^th are the two contributions of astrophysical interest. 
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2.2. T and 4>t{v) calculations 

We now consider how to measure the kinetic temperature T from the complex profiles of a 
pair of elements with different masses. It is a task to deconvolve the thermal and the turbulent 
profiles both of which are unknown a priori. We are mainly interested in cases when line profiles 
are dominated by non thermal broadening, (Xt/^th ^ 1) otherwise the problem is trivial. 

Let us assume that two ions have different masses, m L , and mi > m\. For a fixed value of v 
within the line profile the optical depth is given according to by 

/+oo 
<t> l tiv')^{v-v')dv' , t = l,2, (9) 
-co 

where u L = k^N^. 

This equation is a convolution of two functions. Its Fourier transform is 

/+°° 
r L (v) e 2 ™^ dv = u L <j>i h (q) ,4=1,2, (10) 

-co 

where (j>^(q) an( ^ ^til) are the Fourier transforms of the line profile function <^ h (u) and the 
radial- velocity distribution function 4>[(v), respectively. 

If the chosen ions trace the same volume elements and the ionization fraction is equal for 
both of them (this assumption may not be universally true for any pair of ions but is expected to 
be valid for ions with similar ionization potentials), then <ftl(v) = 4>1(v). In this case two Fourier 
transforms f\{q) and ^(q) yield the exact equations relating the line shapes in the Fourier space 
to the kinetic temperature independently of the velocity and density fluctuations along a given line 
of sight : 



3f? 



Ml) 



0, (11) 



where a = 2/cb(w<2 — mi)/{m\m2) is a numerical constant. 

In the calculation of T from equation ([ll] ) the first step is to determine the Fourier transforms 
ti and T2- The most appropriate way would be to represent r by a set of functions which have 
analytical Fourier transform. We have chosen Gaussian functions since irregular Doppler shifts of 
the local absorption coefficient (^) along the sightline, yielding the apparent complexity of the line 
profiles, result in a mixture of Gaussians (||). The line profiles are modeled using the following 
sum of Gaussians 

, M 

= -7=r £ W iV ~ Vj) /P? , * = 1, 2 . (12) 

V 71 Pi- j= i 

Here M is the number of modes and a\j and a,2,j are their weights. Since basic assumptions of our 
method lead to the similarity of both profiles, we may set the number of modes and their centers 
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to be identical, i.e. M\ = Mi = M and v\j = 1)2 j = Vj. All modes describing the individual profile 
have the same dispersion, denoted fix for the first profile and P2 for the second, with (3% < (3i. 
The choice of the equidispersion Gaussians is by no way restrictive in our case because the modes 
in ( |j~2| ) are used only to describe the entire shape of each line profile and do not bear any physical 
sense as it is in the Voigt profile fitting procedures. The number of modes M is an arbitrary 
parameter and is determined by a satisfactory fit of t\ and T2 solely. The equidispersion modes 



with identical positions allow very easy to meet the second condition in (|llj), - zero imaginary 



part of the ratio f^/fi, - resulting in the requirement axj/a^j = £ = constant for all j. 



It should be noted that the multi-Gaussian model (|12| ) is not the unique one to describe the 
shape of profiles. Any other representations can be used as well. Nevertheless the model ( |i~2|) 
has been favored because of its clearness and easiness in realization. The similar profiles can be 
always fitted by using this model. And on the contrary, if the fitting according to (^) fails then 
the profiles are non-similar and hence our method cannot be applied. 

Equation ( |i~2| ) shows that the total number of free parameters to be found is equal to 2M + 3, 
and that uj can be calculated as 

M 

w t = XX.i ' 1 = 1 ' 2 ' ( 13 ) 

with = 

Using expansion (|l2|), the Fourier transforms of the line profiles can be expressed as 

M 

f L {q) = e-^^i 2 a*J ^ q ,6 = 1,2. (14) 
3=1 



From equations (|i~l|-14) it follows that the kinetic temperature is determined by 

T=((3 2 1 -(3l)/a. (15) 



Having estimated the kinetic temperature, the radial- velocity distribution function 4>t(v) can 
be calculated through the inverse Fourier transform 



of the function 



t(9) 



Mq)e- 2mvc, dq 

-00 



1 = 1,2 



wi <t> l th (q) 

From equations ([14]), (Il6t) and (|T7|) we can readily show that 



H M 

= / ,*2 2 N ^'J eX P 



(v- Vj ) 2 



i = 1 or 2 



(16) 
(17) 

(18) 
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Although in our model f3± and @2 are arbitrary parameters without any physical meaning, their 
natural lower limits should not be smaller than the corresponding thermal widths, otherwise one 
may be trapped into the fitting of noise. Then it follows that the ratio £ = P2/P1 ranges between 
\Jm\/m2 and 1 yielding 01 — v 2 h L > for any parameter sets. 

3. Practical implementation 

In the foregoing section we formulated a method to derive the kinetic temperature by using 
the Fourier transform technique. We have so far dealt with noiseless data. Now we turn our 
attention to the case of noisy data. 

To calculate T from equation (|l^) we need first of all to estimate 2M + 3 parameters of 
model defined by equation ( |l2|) : pi, C> £> v ji °lj> where j = 1,...,M. Since we are primarily 
interested in the recovery of T, another set of the parameters appears to be more convenient : 
9 = {T, C, £, Vj, aij} = {81, ...,#2M+3}> with standing for a parameter vector. The standard x 2 
minimization is used to estimate 9 : 

i=l 

where is the simulated (or observed) normalized intensity within the ith pixel of the line 

profile and I*(Xi,8) = exp{— r(Aj, 0)} is the calculated intensity. Both are convolved with the 
instrumental function <j) sp . Further, A/j is a simulated (or experimental) error, M' is the number 
of pixels involved in the sum, and v = M' — 2M — 3 is the number of degrees of freedom. As 
already mentioned above the choice of the number of modes M is somewhat arbitrary but the 
minimum number of M is determined by the appropriate quality of the fit (x 2 <^ 1)- 

Model calculations with different synthetic spectra have shown, however, that the minimization 
according to (|l^) does not allow to recover the kinetic temperature with a sufficiently high 
accuracy (say 10%). The deviations of the recovered temperature from its underlying value 
exceeded in some cases 100%. This failure is easy to understand : our model for I* is constructed 
using several smoothing operators [convolution with the instrumental function and summing in 
(|l~2|) 1 and, as a result, we have here a typical ill-posed problem (see, e.g., Tikhonov Sz Arsenin 
1977). Parameters in ill-posed problems are known to be very sensitive to data perturbations, 
which results in the inappropriate large dispersions. The ill-posed problems often occur in the 
astronomical image reconstructions as well, and several methods have been developed to overcome 
this obstacle (see Wu 1997 and references cited therein). In order to stabilize the estimation of 
9, we need to augment % 2 by a regularization term penalizing large values of parameters' errors. 
Applied to our particular problem this technique is described below. 



AL 



(19) 
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3.1. The Entropy-Regularized x 2 -Minimization procedure 



This section concerns the Entropy-Regularized x 2 -Minimization [ERM] procedure developed 
to estimate the best fit value for the parameter vector 6. To begin with, we formulate the ERM 
method. Then we state the problem in mathematical language. Finally, we discuss some practical 
implementations. 



where ip is a penalty function and a the regularization parameter which is to be estimated. By 



The choice of ip, especially for nonlinear problems, occurs rather heuristically. Since our 
purpose is to get the most accurate estimation of T, it is reasonable to make ip penalize large error 
in temperature. Having aT = 0^(1 — £ 2 ) and AT/T oc l/fti (1 — C 2 ) t ne easiest choice will be 



which means that a is in units of velocity. 

In spite of its simplicity this function works quite effectively and enables us to locate the true 
value of temperature very closely. [It should be noted that when we change the parameter set 
and make T be a free variable, the form of the regularization term remains the same but with 
f}\ = y/ aT '/(l — Q 2 ). This corresponds to the minimization of the relative error in f3±\. 

The determination of the optimal value of a is the most crucial point in all regularization 
procedures. It remains an open question still for linear problems, although the method of Lagrange 
multipliers, Bayesian estimations and maximum entropy (MaxEnt) concept proved to be very 
successful in many cases. Optimality criterion applied for our problem is motivated by the idea of 
minimization of the cross-entropy which is in principle equivalent to the MaxEnt approach (Wu 
1997). In detail our procedure is described as follows. 

Since entropy is a probabilistic measure of the uncertainty of the information obtained, we 
will consider the parameter vector 6 as a random variable with the probability density function 
(pdf) p a (0). In this notation index a marks the pdf which corresponds to a given value of the 
regularization parameter a. Although p a (0) may have in general a skew peak, we will assume that 
a Gaussian distribution centered at 9q captures most of the probability mass of the true pdf. Then 




minimizing (^) we find the most probable value of the parameter vector 9q. 



1 





1 



e W {--(6-9 )C- 




^det(2^C) 



where C is the error-covariance matrix. 
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This matrix may be estimated through the matrix of second derivatives of C a (8) according to 
its expansion into a Taylor series in the vicinity of the minimum : 

d 2 C a (6 c 



d9 2 



AO 1 



and 



ld 2 C a ( 



2 ddiddi 



(23) 
(24) 



If the errors are uncorrelated, equation ( [22] ) simplifies to the product of the individual 
parameter distributions : 



Pa(9) 



1 V 



(25) 



i=i 



where v' = 2M + 3 and Oj are the parameter dispersions. 



To estimate exactly the level of C a that corresponds to the region containing a certain 
percentage (say 68% = la for normal distribution) of the total probability we would need 
to carry out quite extensive Monte-Carlo simulations at each step of our iterative procedure. 
However, in order to estimate roughly the parameter dispersions we may use instead of 
AC a = £ a {9o + A0) — C a {0o) the corresponding value of A% 2 / which is computed according to 
Press et al. (1992). Since the regularization term in C a is kept small as compared with the x 2 
term, this estimation of the dispersions is expected not to deviate significantly from the true value. 
Then the variances in (|25|) can be evaluated as 



a ; 



(26) 



The definition of the cross-entropy which is the distance between the two distributions p± and 
P2 (called also Kullback measure) is given by (cf., Csiszar, 1996; Wu, 1997) : 

(27) 



fc(pi\\P2) = / Pi (a?) log[pi(x)/p 2 (x)] dx > 0. 
It is the non-additive, non-symmetric measure being equal to only in case p\ 



P2- 



Now suppose that the optimal pdf p O pt(0) is known. Then to estimate the appropriate a- value 
we have to minimize the distance between p opt and p a , i.e. the cross-entropy K.(p a \\p pt) • Taking 



the pdf in the form (25), we obtain 

0- 2 (a op t) 



£(Pa||P. 



opt; 



a] (a) 



1 

+ 2 



o- 2 (a op t) 



+ 



[Oojja) - Oo,j(a op t)Y 
2a](a opt ) 



(28) 



In general we are not given a priori the optimal distribution, but we can construct it adaptively 
from the sequence of solutions obtained for successive values of a ranging from a m i n to a max with 
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a step of Aa. Using the estimation from the previous step as a current approximation for p opt we 
obtain the following iterative sequence : IC(p ai \\p amin ), K,(p a2 \\p ai ) , ... ,lC(p ainax \\p ainax _ Aa ). 

Numerical experiments with synthetic spectra have shown that this sequence really converges 
to 0. The optimum value of a op t, i-e. the value for which the recovered T* was most close to 
the underlying T used to produce the spectrum was found to coincide with the point where the 
curvature of the trajectory /C(oj) = /C(p Qi | I^q,^) passed through the maximum (cf. Lucy 1994). 
Fig. 1 illustrates this behavior. 

As known, the curvature of the continuous curve y{x) is 

ix = y/(ij 2 + l) 3 / 2 , (29) 

where dots indicate differentiation. It is obvious that the value of fi depends on the units in which 
y and x are measured. Therefore in order to get uniform results we always computed \x for the 
values JC and a normalized to 1, i.e. in the coordinates 

K,,Q. = (K, - /Cmin)/ (/C max - /Cmin) , (a - a m in)/ (amax - Omin) • (30) 

The derivatives in (|2^) are computed numerically through the common finite differences: 
fCi = (/Cj + i — 2/Cj + /Cj_i)/Aa 2 and /Q = (/Q+i — /Cj_i)/2Aa. Since numerical differentiation 
is rather sensitive to any even small uncertainties in the function values, to get the smooth 
distribution of [i we need to calculate the K, values with very high accuracy. This means that 
the minimization of £ is to be continued till the stable global minimum with |A£/£| < 10 -6 is 
achieved. The approximation of KL(a) with smoothing spline can be used as well to obtain the 
stable numerical derivatives. But it should be noted that if /C(a) is computed with significant 
errors, then the smoothing spline may noticeably move the maximal curvature point to a wrong 
value of a op t- Thus the only way to derive the reliable estimation of a op t (and hence of T) is to 
carry out the calculations on each step of the procedure as accurate as possible. 

It should be pointed out before going ahead that we can test the estimated value of a op t 
by using the intrinsic properties of the Kullback measure. As noted above, cross-entropy is 
the non-symmetric measure, that is, K,{p\\\p2) 7^ K,{p2\\pi). This property can be used to 
verify the obtained estimations in the following way : we compute /C-values for both upward 
a = Omin (Aa) a max and downward a = a max (— Aa) a m \ n sequences and then estimate the 
corresponding optimal values of a. Their coincidence for both directions may be considered as a 
simple validity proof. Another useful property is the non-additivity of the Kullback measure, i.e. 
fc{Pa+Aa\\Pa) / ^(Pq+Aq/2 I \Pa) + fc(Pa+Aa | IPa+Aa/2) • It means that if the solution of the inverse 
problem does exist then the /C(a)-sequences should yield the same a opt -value independently of the 
step size Aa. 

The limiting values of a m ; n and a max should be chosen in such a way that the expected 
goodness-of-fit does not deteriorate strongly (x 2 ^1)- Although the natural choice of the minimum 
value seems to be a m \ a = 0, in some cases the a-sequence may be continued into the region of 
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negative a keeping inequality C a > valid. The extension of the a-sequence below zero becomes 
especially essential step when the preliminary calculations reveal a low gradient /C(a)-curve over 
the whole interval [0, a max ] indicating that the maximum curvature point lies outside. 

The method described makes it possible to recover the underlying kinetic temperature from 
spectra of different complexity and quality. It should be mentioned, however, that the procedure 
yields only a 'best fit' value for parameter(s) without confidence intervals (a point-like estimation) 
since the dispersions obtained can be over- or underestimated due to the reasons considered above. 
The actual error limits may be revealed only through the numerical experiments with synthetic 
spectra. The model calculations with noisy data (SNR = 100) and spectral resolutions (FWHM) 
ranging from 3 km s _1 to 12 km s _1 have shown, however, that similar unsaturated profiles of C + , 
Si + and Fe + ions provide the accuracy of the kinetic temperature estimation as high as 10%. 



In the previous section we have considered the ERM algorithm. The remaining problem now 
is to show the accuracy of the ERM solutions and hence their reality. 

Since our computational procedure contains several ad hoc assumptions which can hardly be 
backed with unambiguous mathematical justification, it is natural to prove its validity through 
the recovering the kinetic temperature from the synthetic absorption-line spectra at first. For this 
purpose we have generated a variety of line profiles employing different sets of physical parameters 
and several random realizations of the one-dimensional velocity fields. The calculations were 
performed for the pairs of lines C IIA1334 + Si IIA1260 and C IIA1334 + Fe IIA2600 which are 
frequently observed in interstellar and quasar absorption spectra. The procedure may be applied 
to any other suitable pairs as well. 

The physical parameters specifying our cloud model are the column densities of the elements 
in a pair, the dispersion of the turbulent velocity field at, and the kinetic temperature T. Besides, 
all synthetic spectra have been convolved with the instrumental profile of different FWHM. In 
order to make the synthetic spectra look like true observations with an SNR of 100, corresponding 
Gaussian noise has been added. The model parameters are listed in Table 1. 

To calculate synthetic spectra we need to generate stochastic velocity fields at first. Modeling 
of the stationary random velocity fields v(s) with different correlation ranges in space is performed 
through the moving average method described in Appendix A. To generate the correlated fields, 
the following parameterized family of correlation functions was chosen (see Elliott et al. 1997) : 



where t is a linear scale proportional to the distance over which there is appreciable correlation 
between the velocity values at two points (see Monin & Yaglom 1975). This family is widely used 



3.2. Numerical experiments 




(31) 
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in the turbulent flow modeling and provides a diversity of scaling behavior as the parameter e 
varies in the range — oo < e < 2. For instance, regions of negative correlations for the velocity field 
may be obtained with e < 1, whereas for e > 1 we have the regime of positive correlations with 
very long ranges. Parameters i and e used to generate the velocity fields are listed in Table 1 as 
well (£ is given in units of L). 

Listed in Table 1, models 1-11 were constructed to settle the accuracy of the kinetic 
temperature measurement at different spectral resolutions. The synthetic spectra were calculated 
for the fixed values of T = 15000 K and a t = 15 km s _1 , using different velocity fields. The 
velocity fields were chosen rather at random because our aim was mainly to show the ability of the 
method to recover the correct temperature from almost any velocity fields. Thereafter the spectra 
were convolved with Gaussian functions with FWHM = 3, 6, 7, 9, and 12 km s _1 . The spectral 
resolution of 3 km s _1 was achieved in the studies of the interstellar absorption lines with HST, 
whereas the intergalactic spectra are being observed at present with the resolutions 6-12 km s _1 . 
The ratio FWHM/^th2 m our model series lies within the range 1-4 meaning that the smoothing 
effect for FWHM > 6 km s _1 becomes very pronounced. 

To apply our ERM procedure to another kind of absorption spectra, we have calculated the 
synthetic spectra for cr t = 5 km s _1 and the temperatures T = 15000 K and 2000 K (model 12 
and 13, respectively). The last two models 14 and 15 were constructed to test the procedure for 
different ratios of m\jmi and o" t /fth,2- 

Table 2 lists the recovered kinetic temperatures T* along with some other estimated 
parameters. The comparison of columns (7) and (8) shows the high efficiency of the proposed 
computational procedure which allows to estimate T with an accuracy of \AT/T\ <^ 10% even for 
the rather smeared out (models 7-11) spectra. 

In column (6) of Table 2 we list the values of £ = (ih/Pi- Numerical experiments have shown 
that for C ^ 0.9 the dispersion of T*, or*, calculated through ( |26| ) changes its behavior and begins 
to increase with increasing the regularization parameter a. Such peculiar pattern of <tt* stems 
from the fact that at £ ;> 0.9 the off-diagonal elements of the error-covariance matrix C in ( p2[ ) 
are already not small as compared with Cu and, hence, the approximation (25) becomes rather 



vague since the errors become strongly correlated. However, the iterative procedure continues to 
deliver an acceptable T estimation. The procedure fails when £ ~ 1. This is the case when the 
ratio FWHM/w t h,2 is large and the intensity fluctuations within the line profiles are considerably 
smoothed which in turn makes the widths of the subcomponents practically indistinguishable. 

The consequent steps of our computational scheme are illustrated in Figures 2-4. The case 
of positive correlations for the large scale velocity field (e = 1.1 and 1.0) is presented in Fig. 2 
and 4, respectively. In Fig. 2, panel (a) shows with the solid line the adopted correlation function 
Q i£ , and with the dotted one - the autocorrelation function (ACF) calculated for the random 
realization of v(s). The computed ACF may deviate from the exact correlation function because 
of the finite thickness of the absorbing region. Panel (b) shows the corresponding velocity field 
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u(x) = v(x)/at, where x = s/L. This field leads in turn to the complex profiles of C II and Si II. 
Both profiles are convolved with a Gaussian function having the width of 3 km s _1 (FWHM), and 
then a Gaussian noise of S/N = 100 has been added. The final patterns are shown in panels (c) 
and (d) by dots with corresponding error bars. The C II and Si II lines were fitted simultaneously 
employing (|l2|), ( |l9| ) and (|20|). The fit was performed with the increasing number of modes M 
until a x 2 - value of ^ 1 was achieved. In panels (c) and (d) the best-fitting profiles with M = 7 
equidispersion components (located at the velocities shown by the tick marks) are plotted as solid 
curves superimposed on the simulated data. The joint x 2 -value per degree of freedom (y = 129) 
is indicated at the bottom of panel (d). Panel (e) illustrates that if a feasible MaxEnt solution 



does exist, then the Kullback measure lC(a) calculated through (28,30) converges with increasing 
a. The optimal value of a = 0.0 corresponding to the maximum curvature of the trajectory 
fC(a) leads to the estimation T* = 15097 K. Given the kinetic temperature, we may calculate the 
radial- velocity distribution function 4>t(v) using equation ( |I8| ) and either C II or Si II spectrum. 
For both sets of data the solution is expected to be identical. The results obtained are shown 
in panel (f) by solid curves. The histogram plotted in panel (f) is the underlying distribution 
function p(v) which corresponds to the velocity profile v(x) = atu(x) shown in panel (b). Since all 
curves are very much alike, we consider this result as an additional justification of our procedure. 

The regime of negative correlations (s = —0.5) for the random velocity field is illustrated in 
Fig. 3 (here FWHM = 7 km s _1 ). The case for the lowest resolution of 12 km s -1 is demonstrated 
in Fig. 4. The fitting parameters for these examples are listed in Table 2. Analyzing all results, 
we may conclude that the procedure yields stable estimations of T and 4>t(v) for the whole range 
of spectral resolutions considered. 



4. Discussion 

4.1. Comparison with the standard Voigt deconvolution 

Based on the Voigt fitting, most previous approaches extracted information on the gas kinetic 
temperature from the values of b found for the various subcomponents of the complex profiles. 
As a result, absorbing regions have usually been characterized by a variety of T values which, 
in some cases, ranges from K up to « 10000 K (e.g., Spitzer & Fitzpatrick 1995). It seems 
quite natural to apply the standard Voigt fitting procedure to our synthetic spectra to compare 
its result with the estimation obtained by the proposed method. 

To carry out such analysis, we generated the C IIA1334 and Fe IIA2600 absorption profiles 
(Fig. 5, model 14 from Table 1) convolved with the instrumental function of FWHM = 3 km s _1 . 
The corresponding Gaussian noise with an S/N of 100 was added as well. Then we considered 
the spectra shown in Fig. 5 as 'observational' data and tried to analyze them in the standard 
manner, i.e. to fit them with Voigt profiles until a satisfactory result with x 2 ^ 1 was achieved. 
A simultaneous fit to the 10 subcomponents gave the reduced \ 2 P er degree of freedom of 0.687 
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{y = 172). 

The parameter estimates are listed in Table 3 along with the corresponding values of the 
gas temperatures T* and the turbulent velocities a* for the 'individual cloudlets'. The revealed 
values clearly demonstrate a wide spread of the apparent 6-parameters resulting in the large 
variations of T : from ~ 10 K to ~ 18000 K whereas the spectra have been produced by using a 
single underlying temperature T = 10000 K. Moreover, taken at face values, the relative widths 
of bi(C II) and 6«(Fe II) lead to unlikely small values of cr* « 1 — 4 km s _1 , i.e. to a significantly 
lower turbulent velocity than adopted a t = 20 km s -1 . 

The same synthetic spectra having been treated by the proposed procedure yield a single value 
for the gas temperature of 10461 K (see Table 2). The corresponding equidispersion deconvolution 
(with M = 11 components) is plotted in panels c and d in Fig. 5. 

These results lead to the conclusion that the spread of the gas temperature found in some 
regions of the ISM and in the extragalactic intervening clouds may be unphysical and stems purely 
from an inappropriate method of data analysis. 

4.2. The C II and Si II profiles observed in QSO 1937 1009 

We turn now to a brief discussion of the application of the present computational technique to 
recover the kinetic temperature from observational data. For this purpose we have chosen the high 
quality Keck spectra of the quasar Q 1937-1009 (Tytler et al. 1996) showing metal absorption 
lines at the redshift z = 3.572. Since this system is the only one where the deuterium absorption 
was observed at high redshift in both the D Lya and Ly/3 lines, it has been thoroughly studied 
using different methods (Tytler & Buries 1997; Levshakov et al. 1998a; Buries & Tytler 1998). 

According to Tytler & Buries, all metal absorption lines seen at this redshift (C II, C IV, Si II, 
Si IV) 'appear to have the same profiles'. This may lead to the suggestion that these metals are 
similarly distributed in the velocity space and hence a single gas temperature could be sufficient to 
describe their profiles. To estimate its value, we have selected the C IIA1334 and Si IIA1260 lines 
observed with FWHM = 8 km s _1 and S/N = 56 and 100, respectively. The original normalized 
data are shown in Fig. 6 by dots with the corresponding error bars. 

The Voigt fitting analysis previously performed by Tytler et al. (1996) revealed two gas 
clouds separated by 15 km s _1 (the blue and the red components). The physical parameters 
found for these components are the following (Tytler & Buries) : Tbi ue = 17100 ± 1300 K, 
Obiue = 3.35 ± 0.56 km s -1 , T rcd = 19000 ± 5800 K, cr red = 6.75 ± 0.28 km s" 1 (where 

0- = frturb/v^)- 

This system has also been analyzed by different method based on the mesoturbulent concept 
(Levshakov et al. 1998a). In this study the H+D Lya, Ly/3 and higher order Lyman series 
absorption lines (up to Ly-19) have been treated using the reverse Monte-Carlo (RMC) procedure 
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(computational details are given in Levshakov, Kegel & Takahara 1998a, b). The gas temperature, 
T*, and the ensemble turbulent velocity dispersion, cf*, estimated within the framework of 
the one-cloud model were found to lie in the ranges from 10900 to 14300 K and from 18.5 to 
22.3 km s -1 , respectively. 

Applying the present procedure to the C II and Si II lines, we obtained the gas temperature 
T* = 14090 K. Table 4 and Fig. 6 illustrate our results. The best-fitting profiles with M = 5 
equidispersion components are shown in panels a and b by solid lines with the tick marks indicating 
the velocity positions for the separate components. Fig. 6c demonstrates the radial- velocity 
distribution function (solid line curve) calculated through equation (|l~8l). The square root of the 
central second moment of this distribution a* = 8.6 km s _1 depicted in panel c may be considered 
as a measure of the one-dimensional velocity scatter along the space coordinate. This quantity 
naturally falls below the ensemble velocity dispersion a* because the sample of the velocity 
components parallel to the line of sight is not complete in the statistical sense (see Levshakov 
<fe Kegel 1997, for details). Also shown in Fig. 6c is a histogram which represents the velocity 
distribution derived from the analysis of the hydrogen and deuterium lines by the RMC procedure. 

In general we would expect H I and metals to give similar temperatures. They would be 
identical if absorbing gas was fully homogeneous. Analyzing our results on the z = 3.572 system 
we may conclude that both RMC and ERM methods give the consistent estimate of the kinetic 
temperature, T* ~ 14000 K. As for the velocity distribution of neutral hydrogen and deuterium 
and the velocity distribution of ionized metals, they apparently differ (see Fig. 6c). This difference, 
however, stems from the observed asymmetry of the deuterium and the optically thin metal 
absorption lines (Tytler & Buries) : the former exhibits the blue-side asymmetry, whereas the 
latter show the red-side skewness. The interpretation of this result may be as follows. It is 
generally believed that most metal absorption-line systems observed in QSO spectra arise in the 
far outer parts of intervening galaxies or in the intergalactic medium where the main source of 
ionization is photoionization due to the background UV-flux. If the temperature of about 14000 K 
is the mean value for the absorbing gas at z = 3.572 which is in thermal equilibrium, then 
variations in the fractional ionization of different species along the line of sight are mainly governed 
by the density fluctuations n gas within this region because the thermal and ionization state of the 
gas is described by the ionization parameter U which is the ratio of the density of ionizing photons 
to the total gas density. Ionization models of optically thin in the Lyman continuum intergalactic 
gas at constant pressure, photoionized by QSOs, show that the equilibrium temperature is weakly 
dependent on U for metallicities below solar and —4.7 < log?7 < —1.8, where T oc U 1 ^ 8 (Donahue 
& Shull 1991). The deviations between the apparent radial- velocity distribution functions may 
occur when the fractional ionization of different elements responds in different ways to the changes 
in U. Such deviations may be utilized to place constraints on the density fluctuations within the 
absorbing gas. This problem will be studied in detail elsewhere. 
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5. Conclusions 

The outcome of our study, aiming to provide a measurement of the gas kinetic temperature 
from similar complex absorption-line spectra of ions with different atomic weights is the following. 

(1) A new method is developed for using interstellar and intergalactic absorption-line spectra to 
investigate the intervening gas temperature. It is based on the assumption that a pair of elements 
has the same temperature and bulk motion within the absorbing region resulting in the similar 
observed spectra. The method utilizes the Fourier transform of the observed profiles with the 
subsequent Entropy-Regularized Minimization of the \ 2 function augmented by a penalty function 
to gain a stable solution of the inverse problem of spectral analysis [the ERM method] . 

(2) The computational procedure has been tested on a variety of synthetic spectra being generated 
in stochastic velocity fields with different correlation scales. The procedure proved to be very 
effective and robust allowing to recover the temperature with high accuracy, \AT/T\ <; 10 %, in 
the range of spectral resolutions from 3 km s~ 4 to 12 km s _1 (FWHM). 

(3) The application of the standard Voigt fitting analysis to the synthetic spectra calculated in 
the framework of the homogeneous temperature model with T = 10 4 K revealed a wide scatter 
of the gas temperature among the individual 'cloudlets', - from 10 K to 1.8 x 10 4 K. This result 
demonstrates that the scattering of T reported in the literature may have non-physical nature and 
may be caused by an inappropriate methodology. 

(4) The asymmetric C IIA1334 and Si IIA1260 lines observed with the Keck I telescope at 
the redshift z = 3.572 towards quasar 1937-1009 by Tytler et al. have been treated by the 
proposed method in order to evaluate the underlying gas temperature. The obtained estimation of 
T = 14090 K lies in the range consistent with the previously derived values of T from the analysis 
of the H+D Lya and Ly/3 and higher order Lyman series profiles by the reverse Monte-Carlo 
method (Levshakov et al. 1998a). 

It should be emphasized once again that the proposed method can be applied to the similar 
spectra only. In reality different (non-similar) profiles are observed as well (see, e.g., Huang et 
al. 1995, or Songaila 1998). They may be caused by a variety of physical reasons, - overlapping 
of cold and hot gas clouds due to the Doppler shifts of their radial velocities, variations in the 
ionizing flux produced by local sources and/or collisional ionizations in higher density sub-regions 
both leading to the diversity in the ionization fractions for different elements, etc., - and hence 
there cannot be a single physical model to treat such observational data. 

The authors are grateful to David Tytler for sharing his calibrated Keck I/HIRES echelle 
spectra of Q1937-1009 and for a careful reading the manuscript. We also thank Wilhelm Kegel 
for his valuable comments and suggestions which help us much in this study. S.A.L. gratefully 
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A. The moving average method for simulating Gaussian random fields 

Gaussian random fields arise in a wide variety of astrophysical problems including gas 
flows, large scale motions, turbulence etc. To study these problems one needs to model the 
stochastic fields with given characteristics (power spectrum or covariance function). The moving 
average method employed in our modeling the intervening velocity fields is able to reproduce the 
true covariances over the larger scales as compared with the most common spectral (Fourier) 
approximation (Pen 1997; Elliott et al. 1997). The description below follows in general that 
presented in Elliott et al.. 

Consider the Gaussian random field v(s) with the covariance function R(s) and the power 
spectrum S(k). R(s) and S(k) are related through common equations representing the direct and 
inverse Fourier transforms correspondingly : 

/+0O 
R(s) e 27riks ds , (Al) 
-oo 

and 

/+oo 
S(k) e~ 27Tiks dk . (A2) 
-oo 

Define now the function G(s) which is the Fourier transform of the square root of the power 
spectrum S(k) : 

/+oo 1 
S 2(k)e 2niks dk . (A3) 
-oo 

Then v(s) can be represented in the form : 

/+oo 
G(s' - s) dW(s') , (A4) 
-oo 

where W{s) is the Wiener process with unit variance (see, e.g., Cox & Miller 1965). 

The moving average method stems from the discretization of (|A4|), resulting in the 
approximation 

J2 G( n -s)z(n), (A5) 



v 

v(s) ^_ 

i=—p 

r2\ 



where Z is a Gaussian random variable with mean zero and variance (Z ) = Ar, and Ar being 
the step of the equally spaced grid points. 

Gaussian fields are considered to be equal if they have the same mean and covariance function. 
Therefore to test whether the approximation ( |A5| ) represents the real field, we have to compare 
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the covariance function of the approximated field with that which was adopted. In this way we 
may estimate the appropriate step size Ar and the number of the grid points p. The proper values 
of Ar and p should provide minimal difference between the calculated and adopted covariance 
functions. 
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Fig. 1. — Convergence behavior of the cross-entropy defined by (p8|). The filled circles connected 
by dotted line show the normalized values of K, calculated in accord with (|30|). Also plotted is 
the curvature /u/ju max (filled circles connected by dashed line) of the trajectory fC{a). The optimal 
value of a op t = 0.15 corresponds to the maximum curvature at point a = 0.167. This example is 
calculated for model 3 from Table 1. 

Fig. 2. — (a) Autocorrelation functions : solid curve - adopted model, dotted curve - calculated 
from the random realization of the velocity field shown in panel (6) . (c, d) Velocity plots of 
the simulated spectra (dots and la error bars) and the calculated profiles convolved with the 
instrumental resolution of FWHM = 3 km s" 1 (solid curves), corresponding to model 2 in Table 1. 
(e) Cross-entropy as function of the normalized regularization parameter a (filled circles connected 
by dotted line) and the curvature of its trajectory (filled circles connected by dashed line). (/) 
Simulated one-dimensional velocity distribution (histogram) and its estimation through the C II 
(or Si II) profiles. See text for more details. 

Fig. 3. — As Fig. 2, but for different velocity field and FWHM = 7 km s _1 (model 4 in Table 1). 

Fig. 4. — As Fig. 2, but for different velocity field and FWHM = 12 km s _1 (model 11 in Table 1). 
The original (unconvolved) spectra are shown in panels (c) and (d) by dotted curves as well. Note 
the doublet-like structure at Av = 40 km s^ 1 in the Si II spectrum which is completely smeared 
out after convolution. 

Fig. 5. — (a, b) Velocity plots of the C IIA1334 and Fe IIA2600 synthetic spectra (dots with la 
error bars, model 14 in Table 1). The smooth curves show the Voigt profiles convolved with the 
instrumental resolution of 3 km s _1 (FWHM) which produce the best simultaneous fit to all data. 
The velocity positions for each of ten absorption components are labeled by the tick marks at 
the top of each panel. (c,d) Same as (a,b), but for the ERM result based on the equidispersion 
deconvolution of the C II and Fe II profiles with a single kinetic temperature (see text for more 
details). 

Fig. 6. — (a, b) Normalized data points with la error bars for HIRES echelle spectrograph 
observations (FWHM = 8 km s _1 ) on the Keck telescope of the C IIA1334 and Si IIA1260 absorption 
features (the signal-to-noise per pixel is, respectively, 56 and 100) at the redshift z = 3.572 towards 
the quasar 1937-1009 (Tytler & Buries 1997). The ERM result (reduced \ 2 = 0.331 with 24 
degrees of freedom) for C II and Si II (solid lines) is based on simultaneously fitting the multiple 
equidispersion lines located at the velocities shown by the tick marks at the top of each panel. 
The evaluated kinetic temperature T* = 14090 K is plotted in panel a. (c) The radial- velocity 
distribution function <pt{v) restored through the observed profile of the C II (or Si II) line using 
eq. ( |i~8| ) and the model parameters from Table 4. The square root of the central second moment 
a* of the cj>t(v) distribution is also shown in the panel. For comparison, histogram reproduces the 
velocity distribution function estimated by the reverse Monte-Carlo [RMC] procedure (Levshakov 
et al. 1998a) through the analysis of the H+D Lya, Ly/3 and the higher order Lyman series lines 
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(up to Ly-19) from the same absorption-line system. 
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Table 1. Model parameters used to generate synthetic spectra of C IIA1334 and Si IIA1260 

(models 1-13) and synthetic spectra of C IIA1334 and Fe IIA2600 (models 14,15). The 
parameterized family of velocity correlation functions (|3l"l ) is defined by the parameters £ and e 

(see text). 



model 


pair 


N, 


0t, 


T, 


FWHM, 


I 








cm -2 


km s _1 


K 


km s _1 






1 


C II 


4(13) 


15 


1.5(4) 


3 


0.1 


1.0 




Si II 


4(12) 












2 


C II 


4(13) 


15 


1.5(4) 


3 


0.1 


1.1 




Si II 


4(12) 












3 


C II 


4(13) 


15 


1.5(4) 


6 


0.1 


1.1 




Si II 


4(12) 












4 


C II 


4(13) 


15 


1.5(4) 


7 


0.25 


-0.5 




Si II 


4(12) 












5 


C II 


4(13) 


15 


1.5(4) 


7 


0.25 


0.5 




Si II 


4(12) 












6 


C II 


A( 1 X\ 


15 


1.5(4) 


7 


0.25 


— 1.5 




Si II 


4(12) 












7 


C II 


4(13) 


15 


1.5(4) 


9 


0.1 


1.0 




Si II 


4(12) 












8 


C II 


4(13) 


15 


1.5(4) 


9 


0.1 


1.1 




Si II 


4(12) 












9 


C II 


4(13) 


15 


1.5(4) 


9 


0.1 


-1.0 




Si II 


4(12) 












10 


C II 


4(13) 


15 


1.5(4) 


12 


0.1 


1.1 




Si II 


4(12) 












11 


C II 


4(13) 


15 


1.5(4) 


12 


0.1 


1.0 




Si II 


4(12) 












12 


C II 


4(13) 


5 


1.5(4) 


3 


0.5 


-1.5 




Si II 


4(12) 












13 


C II 


4(13) 


5 


2.0(3) 


3 


0.05 


-1.5 




Si II 


4(12) 












14 


C II 


5.5(13) 


20 


1.0(4) 


3 


0.3 


-1.5 




Fe II 


5.3(12) 












15 


C II 


5.5(13) 


20 


1.0(4) 


6 


0.3 


-1.5 




Fe II 


5.3(12) 













Table 2. Model parameters derived from the synthetic spectra through the equidispersion 
deconvolution (12) with the subsequent entropy-regularized minimization of the objective 

function ©. 



model 


G^min 


Oopt 


M 


x 2 


c 


T*f K 


T, K 


1 


0.0 


0.2 


8 


0.713 


0.86 


14271 


15000 


2 


-0.4 


0.0 


7 


0.674 


0.77 


15097 


15000 


3 


0.0 


0.15 


4 


0.725 


0.89 b 


13775 


15000 


4 


-0.15 


0.0 


4 


0.364 


0.84 


14576 


15000 


5 


0.0 


0.275 


4 


0.797 


0.85 


15401 


15000 


6 


0.0 


0.2 


5 


0.814 


0.86 


15866 


15000 


7 


-0.037 


-0.032 


4 


0.710 


0.93 b 


14530 


15000 


8 


-0.069 


-0.062 


4 


0.409 


0.89 


13537 


15000 


9 


0.0 


0.012 


4 


0.741 


0.95 b 


13910 


15000 


10 


-0.006 


-0.003 


2 


0.408 


0.96 b 


14473 


15000 


11 


-0.0203 


-0.0198 


4 


0.458 


0.93 b 


14632 


15000 


12 


-0.2 


0.5 


2 


0.257 


0.71 


16112 


15000 


13 


0.0 


0.02 


5 


0.752 


0.92 b 


1816 


2000 


14 


0.0 


0.3 


11 


0.763 


0.75 


10461 


10000 


15 


-0.04 


0.02 


8 


0.756 


0.85 


10081 


10000 



a T* - estimated temperature, T - adopted temperature. 
b increasing dispersion <tt* (see text for details). 
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Table 3. The standard Voigt fitting procedure results. Model parameters evaluated 
simultaneously from the C IIA1334 and Fe IIA2600 synthetic spectra (Fig. 5) by fitting ten 
individual subcomponents numbered in the order of increasing relative velocity Av. Reduced 

X 2 (with 172 degrees of freedom) is equal to 0.687. The adopted values for the kinetic 
temperature, T, and the rms turbulent velocity dispersion, a t , are, respectively, equal to 10 4 K 

and 20 km s _1 (model 14 in Table 1). 



No An bi(C ID 


6i(FeII), 


T* 




km s _1 km s _1 


km s _1 


K 


km s" 1 


1 -47.86 3.75 


2.09 


8896 


0.9 


2 -37.79 6.46 


4.90 


16212 


3.1 


3 -27.42 5.998 


5.997 


11 


4.2 


4 -12.68 5.79 


3.64 


18590 


2.0 


5 -5.90 5.06 


2.50 


17821 


0.7 


6 0.52 5.20 


3.76 


11854 


2.3 


7 9.78 6.07 


4.90 


11837 


3.2 


8 23.24 6.00 


5.37 


6564 


3.7 


9 33.93 5.86 


5.63 


2467 


3.9 


10 41.90 4.38 


3.18 


8310 


1.9 


Table 4. The fit based on the entropy-reg 


;ularized minimization. Model parameters 


evaluated simultaneously from the observed C IIA1334 and Si 


IIA1260 spectra (Fig. 


6) by 


fitting five equidispersion components numbered 


in the order of 


increasing velocity position Vj. 


Reduced x 2 (with 24 degrees of freedom) is equal to 0.331. The estimated value for the kinetic 


temperature is T* = 14090 K, and the widths of the components are equal to 




13* {C II) = 5.735 km s" 1 and /3*(Si II) 


= 4.671 km s- 1 


(see equation (|l2])). 




No. Vj, km s _1 


aj(C II) 


% (Si II) 


1 -15.18 


0.813 




0.833 


2 -7.56 


1.619 




1.658 


3 -0.93 


2.998 




3.071 


4 8.57 


4.930 




5.051 


5 17.62 


0.445 




0.456 
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Figure 2 




0.0 0.2 0.4 -20 20 0.0 0.5 1.0 

scale shift, As/L velocity, km s 1 regularization parameter 




Figure 3 



1.0 




i 1 i 






ACF 






1 = 0.1 


0.5 




6 = 1.0 


0.0 








a 




0.5 




1 , 1 



0.0 0.2 0.4 



scale shift, As/L 




0.0 0.5 1.0 

coordinate, x = s/L 

Figure 4 




velocity, km s 1 regularization parameter 




-40 -20 20 40 60 _ 40 -20 20 40 60 

velocity, km s" 1 velocity, km s" 1 



Voigt Deconvolution 




Figure 5 




Figure 6 



